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The linear stability of MHD Taylor-Couette flow of inflnite vertical extension is considered for 
liquid sodium with its small magnetic Prandtl number Pm of order of 10~^. The calculations are 
performed for a container with r) = 0.5 with an axial uniform magnetic field excluding counter 
rotating cylinders. The sign of the constant a in the basic rotation law Q, = a + b/R^ strongly 
influences the presented results. It is negative for resting outer cylinder. The main point here is 
that the subcritical excitation which occurs for large Pm disappears for small Pm (cf. Fig. 
This is the reason that the existence of the magnetorotational instability remained unknown over 
decades. 

For rotating outer cylinder the limiting case a = (i.e. p, = ff) plays an exceptional role. 
The hydrodynamic instability starts to disappear while the hydromagnetic instability exists with 
minimal Reynolds numbers at certain Hartmann numbers of the magnetic field. These Reynolds 
numbers exactly scale with Pm"'^''^ resulting in moderate values of order 10"* for Pm=10~^. However, 
already for the smallest positive value of a the Reynolds numbers start to scale as 1/Pm leading 
to much higher values of order lO'' for Pm=10~^. Hence, for outer cylinders rotating faster than 
the limit a = it is exclusively the magnetic Reynolds number Rm which directs the excitation 
of the instability. They are resulting as lower for insulating walls ('vacuum') than for conducting 
walls. Generally, the magnetic Reynolds numbers for liquid sodium have to exceed values of order 
10 leading to frequencies of about 20 Hz for the rotation of the inner cylinder if containers with 
(say) 10 cm radius are considered. The required magnetic fields are about 1000 Gauss. 

Also nonaxisymmetric modes have been considered. With vacuum boundary conditions their 
excitation is always more difflcult than the excitation of axisymmetric modes; we never observed a 
crossover of the lines of marginal stability. For conducting walls, however, such crossovers exist for 
both resting and rotating outer cylinders, and this might be essential for future dynamo experiments. 
In this case, however, the instability also can onset in form of oscillating axisymmetric patterns of 
flow and field and the Reynolds numbers of these solutions are lower than the Reynolds numbers 
for the stationary solutions. 

PACS numbers: 47.20.Ft, 47.20.-k, 47.65.+a 



I. INTRODUCTION 

The longstanding problem of the generation of turbu- 
lence in various hydrodynamically stable situations has 
found a solution in recent years with the MHD shear 
flow instability, also called magnetorotational instabil- 
ity (MRI), in which the presence of a magnetic field 
has a destabilizing effect on a differentially rotating flow 
with the angular velocity decreasing outwards. The MRI 
has been formulated decades ago |l|, ^ for ideal Taylor- 
Couette flow, but its importance as the source of tur- 
bulence in accretion discs with differential (Keplerian) 
rotation was first recognized by Balbus and Hawley, 

However, the MRI has never been observed in the lab- 
oratory [@, |[ |[ 0. After Goodman and Ji [|j was the 
absence of MRI due to the small magnetic Prandtl num- 
ber approximation used in [Q. The magnetic Prandtl 
number Pm is really very small under laboratory condi- 
tions (~ 10~^ and smaller, see Table |). 

A proper understanding of this phenomenon is very im- 



TABLE I: Parameters of the fiuids suitable for MHD experi- 
ments taken from |2| and [b| 





P [g/cm^] 


V [cm^/s] 


•q [cm^/s] 


Pm 


Mercury 


5.4 


1.1-10"^ 


7600 


1.4-10"^ 


Gallium 


6.0 


3.2-10"^ 


2060 


1.5-10"^ 


Sodium 


0.92 


7.1-10"^ 


810 


0.88-10"^ 



portant for possible future experiments, including Taylor- 
Couette flow dynamo experiments. The simple model of 
uniform density fluid contained between two vertically- 
infinite rotating cylinders is used with constant magnetic 
field parallel to the rotation axis. For viscous flows the 
most general form of the rotation law Vl{R) in the fluid 
is 

VL{R)^a + b/R^, (1) 
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where a and h are two constants related to the angular 
velocities fiin and fiout with which the inner and the outer 
cylinders are rotating and R is the distance from the 
rotation axis. If i?in and i?out (-Rout > ^in) are the radii 
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FIG. 1: Cylinder geometry of the Taylor-Couette flow with 
axial magnetic field. 

of the two cylinders then 

a=T^^in and 6 = i?f„l^a„ (2) 

1 — 77^ 1 — r/^ 

with 

A = ^out/f^in and fl = Rin/Rout- (3) 

Following the Rayleigh stability criterion, d{R^Vl)^ / dR > 

0, rotation laws are hydrodynamically stable for a > 0, 

1. e. jx> if . They should in particular be stable for rest- 
ing inner cylinder, i.e. /t — > 00. Richard and Zahn fo- 
cused attention to the experimental results of Wendt 
who found nonlinear instability for this case for Reynolds 
numbers of order 10^. The finite-amplitude instability of 
hydrodynamically stable rotation laws must therefore re- 
main in the astrophysical discussion. However, later ex- 
periments with very similar Taylor-Couette flow experi- 
ments for resting inner cylinder demonstrated the results 
of Wendt as due to rather imperfect container construc- 
tions and the flow remained laminar even for Reynolds 
numbers up to 10^, 

One of the targets in the present paper is the axisym- 
metry of the excited modes. We have shown in ||l3| that 
for containers with conducting boundaries it happens for 
sufficiently strong magnetic fields that the mode with the 
lowest eigenvalue (i.e. the lowest Reynolds number) is a 
nonaxisymmetric mode. As an impressive example, in 
Fig. ^ for Pm=0.01 the crossover of the instability lines 
for axisymmetric and nonaxisymmetric modes is shown 
for Hartmann numbers of about 400 (see |20|). Despite of 
its general meaning this behavior is only known so far for 
conducting walls and for magnetic Prandtl numbers not 
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FIG. 2: Instability lines for axisymmetric (solid) and nonax- 
isymmetric modes (m=I, dashed line) for conducting walls 
and Pm=O.OI (f) = 0.5). 

smaller than 10^^ (see |l^). For possible laboratory ex- 
periments we have to extend, however, the computations 
to insulating boundaries (vacuum) and to much smaller 
magnetic Prandtl numbers Pm. 

The equations, therefore, are solved here mainly for the 
single small magnetic Prandtl number Pm = 10~^ very 
close to the value for liquid sodium (see Table ||). The 
aspect ratio of the container walls radii in the present 
paper is fixed to ^ = 0.5. 

II. BASIC EQUATIONS 

The MHD equations which have to be solved are 
dii 1 

— + (uV)u = — Vp + i^Ait + J X B (4) 
at p 

and 

dB 

^— = curl(tt X B)+ rjAB, (5) 
ot 

with the electric current J = cui\B/fiQ and divit = 
divB = 0. They are considered in cylindrical geometry 
with R, (j), and z as the cylindrical coordinates. A viscous 
electrically-conducting incompressible fluid between two 
rotating infinite cylinders in the presence of a uniform 
magnetic field parallel to the rotation axis leads to the 
basic solution Ur — Uz = Br = ^ 0,Bz ^ Bq ~ 
const., and — aR + b/R, with U as the flow and B 
as the magnetic field. We are interested in the stability 
of this solution. The perturbed state of the flow may be 
described by m'^, u'^, u[, B'^, B'^, B'^, p' with p' as the 
pressure perturbation. 

Here only the linear stability problem is considered. 
By analyzing the disturbances into normal modes the 
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solutions of the linearized hydromagnetic equations are 
of the form 

B' = B{R) ei('"<^+'=^-'^*), 

u' = u{R) e'('"'^+*=^-"*). (6) 

From hereon all dashes have been omitted from the no- 
tations of fluctuating quantities. Only marginal stability 
will be considered hence the imaginary part of uj, i.e. 
T{lo), always vanishes. We use 



H — \/ i?in(^out — -Rin) 



(7) 



as the unit of length, the r\lli as the unit of perturbed 
velocity and Bq as the unit of perturbed magnetic field 
and work with the magnetic Prandtl number 



Pm 



(8) 



with V the kinematic viscosity and r] the magnetic dif- 
fusivity. Note also as the unit of wave numbers 

and v/H^ as the unit of frequencies. After elimination 
of both pressure fluctuations and the fluctuations of the 
vertical magnetic field, B'^, the linearized equations are 



duR UR im 



(9) 
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(11) 
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(13) 



Here the Reynolds number Re and the Hartmann number 
Ha are defined as 



Re 



f^inJ^iii(^out — Rin) 



and 



Ha = B, 



01 



' -Rin (-Rout — Rii 



(14) 



(15) 



For given Hartmann number and magnetic Prandtl num- 
ber in the present paper we shall compute with a linear 
theory the critical Reynolds number of the rotation of 
the inner cylinder, also for various mode numbers m. 



III. BOUNDARY CONDITIONS, NUMERICS 

An appropriate set of ten boundary conditions is 
needed to solve the system (^. . .([l^). Always no-slip 
conditions for the velocity on the walls are used, i.e. 
Ur, = U0 = duR/dR = 0. The boundary conditions for 
the magnetic field depend on the electrical properties of 
the walls. The tangential currents and the radial com- 
ponent of the magnetic field vanish on conducting walls 
hence dB^/dR + B^/R = Br = 0. These boundary con- 
ditions may also hold both for R — Ri^ and for R — i?out ■ 

The homogeneous set of equations (||) . . . (|l^ together 
with the boundary conditions determine the eigenvalue 
problem of the form £(fc, m. Re, Ha, 7?.(a;)) — for given 
Pm. The real part of cu, i.e. TZ{uj), describes a drift 
of the pattern along the azimuth which only exists for 
nonaxisymmetric flows. For axisymmctric flows (m = 
0) the real part of w, i.e. TZ{uj), is zero for stationary 
patterns of flow and field and it is nonzero for oscillating 
solutions, which are called overstability. £ is a complex 
quantity, both its real part and its imaginary part must 
vanish for the critical Reynolds number. The latter is 
minimized by choice of the wave number k. TZ^u) is the 
second quantity which is fixed by the eigen equation. 

The system is approximated by finite differences with 
typically 200 radial grid points. The resulting determi- 
nant, C, takes the value zero if and only if the values Re 
are the eigenvalues. We can also stress that the results 
are numerically robust as an increase of the number of 
grid points does not change the results remarkably. For 
a fixed Hartmann number, a fixed Prandtl number and a 
given vertical wave number k we find the eigenvalues of 
the equation system. They are always minimal for a cer- 
tain wave number which by itself defines the marginally 
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unstable mode. The corresponding eigenvalue is the de- 
sired Reynolds number. 

The situation changes for insulating walls. The mag- 
netic field must match the external magnetic field for 
vacuum. It is known for this case that the boundary 
conditions for axisymmetric solutions strongly differ from 
those for nonaxisymmetric solutions (see [p5[). The con- 
dition curl/jB = in vacuum immediately provides 



kR 



(16) 



at i? = i?in and R = i?out- From the solution of the 
potential equation — one finds 



B 



iB, 



(- 



I,n{kR) \kR 
for R ~ i?in and 
\B 



Im{kR)+I,n+iikR)) =0 (17) 



B 



R 



(- 



K,nikR) \kR 



K^(kR) - K,n+i{kR)) -0 (18) 



for R ~ i?out- Im and Km are the modified Bessel func- 
tions (with different behavior at i? — > and R oo). 
One can eliminate with divB=0 the vertical component 
Bz of the magnetic field in the boundary conditions 



IV. RESULTS 
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FIG. 3; Marginal stability lines for axisymmetric modes with 
resting outer cylinder of conducting material. The shaded 
area denotes subcritical excitations of unstable axisymmetric 
modes by the external magnetic field. 



The following results concern different aspects of 
the MHD Taylor-Couette problem for small magnetic 
Prandtl number Pm. In Section A the main realization 
of the case a < (here with resting outer cyhnder, i.e. 
/i = 0) is considered. There is instability even without 
magnetic fields so that the bifurcation lines start at the 
y-axis. In Section B the special case a = is considered 
with very surprising results. The Section C presents the 
results for the two experiments with jl — and ft = 0.33 
with respect to the axisymmetry of the eigenmodes. In 
Section D the existence of oscillating modes is discussed, 
i.e. the case of overstability for small magnetic Prandtl 
numbers. 



A. Subcritical excitation for large Pm (a < 0) 

Figure |^ shows the stability lines for axisymmetric 
modes for containers with conducting walls and with rest- 
ing outer cylinder for fluids of various magnetic Prandtl 
number. Only the vicinity of the classical hydrodynamic 
solution with Re = 68 is shown. There is a strong dif- 
ference of the geometry of the bifurcation lines for Pm 
> 1 and Pm < 1. In the latter case, for fluids with low 
electrical conductivity the magnetic field only suppresses 
the instability so that all the critical Reynolds numbers 
exceed the value 68, and this the more the stronger the 
magnetic field is. 



For sufficiently small magnetic Prandtl number the sta- 
bility lines hardly differ, which is the situation already 
considered by Chandrasekhar without any indication 
of magnetorotational instability. 

The opposite is true for Pm > 1. Note that in Fig. ^ 
for materials with high electrical conductivity the result- 
ing critical Reynolds numbers are smaller than Re = 68. 
The magnetic field with small Hartmann numbers sup- 
port instability patterns rather than to suppress them. 
This effect becomes more effective for increasing Pm but 
it vanishes for stronger magnetic fields. Obviously, the 
MRI only exists for weak magnetic fields and high enough 
electrical conductivity and/or molecular viscosity (when 
the fields can be considered as frozen in and/or enough 
viscosity prevents the action of the Taylor-Proudman the- 
orem) . 

Note that the subcritical excitation of Taylor vortices 
only works for weak magnetic fields. The upper limits 
of the possible Hartmann numbers can be observed for 
the magnetic Prandtl numbers 1 and 10 in Fig. ^. After 
our computations the subcritical excitation of Taylor vor- 
tices for weak magnetic fields requires rather high mag- 
netic Prandtl numbers. The microscopic values for Pm 
are orders of magnitudes smaller than unity, so that there 
should be no chance to realize the subcritical excitation of 
Taylor vortices by experiments. However, the speculation 
may be allowed whether really the microscopic Pm is the 
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FIG. 4: Critical Reynolds number versus i) for ji — rf and 
Pm= 1. 



basic input. The scenario is also interesting whether pos- 
sible finite-amplitude hydrodynamic instabilities provide 
some kind of background turbulence which can be con- 
sidered as modifying the value of the magnetic Prandtl 
number, ||^. The turbulence influences both the viscosity 
values and the magnetic-diffusivity values so that 



Pm — > Pnicff = 



(19) 



with i"T and ryx as the eddy viscosity and the eddy diffu- 
sivity, resp. Because of the existence of the pressure term 
in the momentum equation, both quantities are not iden- 
tical. We do not have precise knowledge about the effec- 
tive turbulent magnetic Prandtl number but it has been 
demonstrated that values of order 0.1 or somewhat larger 
should not be unlikely, jl^. Insofar if such speculations 
are not too far from the reality, it is not completely clear 
that the subcritical excitation of Taylor vortices which 
we have presented in Fig. |^ is unobservable in general. 



B. The case a = (/i = 17^) 

There is a universal scaling on Pm for the special case 
with a = in the basic flow profile ([l]), i.e. for jl = ff' . 
Then the term with d{R'^Vt) /dR in Eq. ( |l0| ) vanishes and 
for m — Lo — Q one finds that the quantities ur^Uz, Br 
and Bz are scaling as Pm^^/^ while u^^B^^k Ha scale 
as Pm°. Then also the Reynolds number for the axisym- 
metric modes scales as 



Re oc Pm 



-1/2 



(20) 



The scaling does not depend on the boundary conditions 
as these for m = also comply with the relations. 

The result ( pO| ) has numerically been found by Willis 
and Barenghi for vacuum boundary conditions, W^- 



However, Riidiger and Shalybkov 
found the much steeper scaling 



Pj for a > (A >rf) 



ReocPm^\ (21) 
resulting in the surprisingly simple relation 



for the magnetic Reynolds number Rm 
riini?in(-Rout - Bin) /rj and 



Ha cx Pm 



'1/2 



resulting in 



Ha* cx const. 



(23) 



(24) 



for Ha* ~ Boy^Ri^{Rout — Bin)/ fJ-oPV^ (Lundquist num- 
ber, see [|l^). In case of small magnetic Prandtl num- 
ber the exact value of the microscopic viscosity is totally 
unimportant for the excitation of the instability. In con- 
sequence, however, the corresponding Reynolds numbers 
for the MRI seem to differ by 2 orders of magnitude, i.e. 
10** and 10^. Insofar, experiments with jl = if' seem to 
look much more promising than experiments with fi > if' . 

Unfortunately, this challenging possibility cannot be 
utilized in experiments. The critical Reynolds number 
for jl = if and Pm=l as a function of 77 is given in Fig. 
^. The total minimum of the Reynolds number is 54.4 for 
?} = 0.27 so that after ( ^ ) one expects the value 1.7-10'' 
for the Reynolds number for Pm= 10~^. Fig. shows 
the behavior of this result in the vicinity of /i = 77 . There 
is a vertical jump from 10* to 10^ in an extremely small 
interval of the abscissa. This sharp transition does not 
exist for Pm=l, it is only due to the very small value 
of Pm. For this case in Fig. ^ the coexistence of both 
hydrodynamic and hydromagnetic instability is also pre- 
sented. The jump profile for Pm= 10~^ in Fig. || (right) 
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FIG. 5: Critical Reynolds numbers for the Taylor-Couette 
flow versus jl for fj = 0.27 and Pm=l (left) and Pm=10~^ 
(right). The curve for the hydrodynamic instability (Ha=0) 
is dashed and the hydromagnetic curve (Ha>0) is solid. The 
dotted line denotes the location of a=0, i.e. ji = if . 



Rm cx const. 



(22) makes it clear that such experiments with jl = if are not 
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TABLE II: Coordinates of the absolute minima in Figs. |6| and 
u\ for rotating outer cylinder {fi — 0.33) 





conducting walls 


insulating walls 


Reynolds number 


2.13 ■ 10** 


1.42 ■ 10'' 


magnetic Reynolds number 


21 


14 


Hartmann number 


1100 


1400 


Lundquist number 


3.47 


4.42 



possible. Even the smallest deviation from the condi- 
tion jji — 1)^ drasticaUy changes the excitation condition. 
For ji smaller than 77^ (negative deviations) the hydro- 
dynamic instabihty sets in and for /i slightly exceeding 
if (positive deviations) the Reynolds number suddenly 
jumps by two orders of magnitudes. 



C. Excitation of nonaxisymmetric modes 





FIG. 6: Insulating walls (vacuum): Stability lines for axisym- 
metric (m = 0, solid lines) and nonaxisymmetric instability 
modes with m = 1 (dashed). LEFT: resting outer cylinder, 
RIGHT: rotating outer cylinder (33%). 

Let us now concentrate to the small magnetic Prandtl 
number for liquid sodium, i.e. Pm — 10^^. We start 
with the results for containers with insulating walls and 
outer cylinders at rest, (Fig. ^). There are then linear 



FIG. 7: Conducting waUs: Stability lines for axisymmetric 
(m = 0, solid lines) and nonaxisymmetric instability modes 
with m = 1 (dashed). LEFT: resting outer cylinder, RIGHT: 
rotating outer cylinder (33 %). 



instabilities even without magnetic fields. For Ha= so- 
lutions for m = (Re= 68) and m — 1 (Re= 75) are 
known, see [|l3|. The axisymmetric mode possesses the 
lowest eigenvalue. This is also true for nonvanishing mag- 
netic field; we do not find any crossover of the instability 
lines for axisymmetric and nonaxisymmetric modes. The 
same is true for containers with rotating outer cylinder 
(Fig. ||b). For growing jl the Reynolds number for the 
hydrodynamic solution moves upwards, reaching infin- 
ity for fi = r)^ — 0.25 (here). The MRI is represented 
by characteristic minima, in our case for p, = 0.33 at 
Hartmann numbers of order 10'^ and Reynolds numbers 
of order 10^. The exact coordinates of the minima are 
given m Table |l[ 

In order to characterize the Hartmann numbers note 
that for liquid sodium 



B = 2.2 



Ha 



-Rout/lOcm 



Gauss. 



(25) 



Hence, for i?out — 22 the magnetic field and the Hart- 
mann number have the same numerical values. With 
ly = 10"^ cm^/s and 17=0.5 it follows from (|l^) and (|l^) 

/in = 64 ^^^^^^^^ Hz (26) 



(i?out/10cm)^ 



for the frequency of the inner cylinder. Hence, a con- 
tainer with an outer radius of 22 cm (see above) and an 
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inner radius of 11 cm filled with liquid sodium and em- 
bedded in vacuum requires a rotation of about 19 Hz in 
order to find the MRI. Following ( p5| ) the required mag- 
netic field is about 1400 Gauss. 
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FIG. 8: The same as in Fig. Ma but with the inclusion of 
oscillating axisymmetric modes (overstability) appearing here 
for lower Reynolds numbers. 

The results for containers with conducting walls are 
given in Fig. 0. Note that the minimal Reynolds numbers 



given in Fig. |7|b are higher than for insulating cylinder 
walls. The influence of the boundary conditions is not as 
small as expected. The main difference between the two 
sorts of boundary conditions, however, is the existence of 
crossovers of the instability lines for m = and m = 1 in 
case of conducting walls. For both resting and rotating 
outer cylinder Hartmann numbers exist above which the 
nonaxisymmetric mode possesses a lower Reynolds num- 
ber than the axisymmetric mode. We have already shown 
the existence of such crossovers for conducting walls for 
1 < Pm < 0.01 in |l^. It is now clear that the occurrence 
of nonaxisymmetric solutions as the preferred modes is a 
rather general phenomenon for containers with conduct- 
ing walls which can become important for the design of 
future dynamo experiments (Cowling theorem). 



D. Excitation of oscillating modes 

There are not only stationary patterns of flow and field 
possible but the instability can also onset in form of os- 
cillating solutions. This effect is called overstability. In 
case of rotating convection between two layers heated 
from below the onset of the instability in form of os- 
cillating solutions even possesses the lowest eigenvalues 
for certain Prandtl numbers, We find a very simi- 
lar behavior for the MHD Taylor-Couette flow between 
conducting cylinders for resting outer cylinder (see Fig. 

. It is a pair of waves traveling in positive and negative 
^-direction. Note that the cylinder considered here has 
no bound in vertical direction. If the cylinder is finite, 
however, the possibility exists that traveling waves might 
be combined to standing waves. 
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